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Abstract. - We lay the theoretical basis for one-bead microrheology with rotating particles, 
i.e, a method where colloids are used to probe the mechanical properties of viscoelastic media. 
Based on a two-fluid model, we calculate the compliance and discuss it for two cases. We first 
assume that the elastic and fluid component exhibit both stick boundary conditions at the 
particle surface. Then, the compliance fulfills a generalized Stokes law with a complex shear 
modulus whose validity is only limited by inertial effects, in contrast to translational motion. 
Secondly, we find that the validity of the Stokes regime is reduced when the elastic network is 
not coupled to the particle. 



Introduction. - In recent years the experimental method of microrheology [1] emerged 
as a powerful tool to monitor the mechanical properties of viscoelastic soft materials [2, 3] 
especially in biological systems [4,5] including cells [6]. The main idea is to disperse micron- 
sized beads into the material and monitor their motion either as reponse to external forces 
(active method) [6,7] or due to Brownian fluctuations (passive method) [5,8]. Whereas in 
the first method, the frequency-dependent response function or compliance is measured di- 
rectly, it is inferred in the second method from the particle's positional fluctuations using the 
fluctuation-dissipation theorem and the Kramers-Kronig relations [5] . 

In this article, we address the particles' rotational degree of freedom and lay the theoretical 
basis to use it in microrheology. To monitor the particles' orientations, they have to be 
anisotropic. And indeed, first experimental implementations of one-bead microrheology with 
either birefringent spherial particles [9] or microdisks [10] do exist. 

The main theoretical problem of microrheology is how the particle compliance measured 
in experiments relates to the viscoelastic properties of the material quantified in its complex 
shear modulus. We will show here, based on a model-viscoelastic medium, the so-called two- 
fluid system [11], that the compliance in the rotational case obeys a generalized Stokes law 
which is only limited by an upper crossover frequency in contrast to the translational motion 
where the validity is restricted to a frequency window [5, 12]. So one advantage of rotating 
particles is that they extend the frequency range that allows a straightforward interpretation 
of the experimental results by the Stokes law. 
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Let US first review the reason for the frequency window in the translational case [5]. As a 
model- viscoelastic medium, we consider an elastic component, e.g., a polymer network, which 
can move relative to a Newtonian fluid. A particle embedded in this medium experiences 
an external oscillating force with amplitude F(lu) and reacts with an oscillating displacement 
described by the amplitude x{uj) — a{uj)F{u!), where a{uj) is the compliance as a function 
of frequency uj. For lu = 0, the displaced particle creates a deformation field in the elastic 
network which includes both shear and compressional contributions. For small lu, the elastic 
network thus has to move relative to the incompressible fluid since the latter only allows 
shear motions. However, the frictional force between both components increases with their 
relative speed. Hence, beyond a crossover frequency cucii elastic network and viscous fluid 
are strongly coupled and can be considered as a single incompressible viscoelastic medium. 
The compliance for a spherical particle then assumes the simple generalized Stokes relation 
a(w) = [67rG(w)a]~^ where G(w) = ^ — iur] is the complex shear modulus, rj is viscosity and /i 
denotes Lame's elastic constant associated with shear. Beyond a second crossover frequency 
ujc2, inertial effects of the fluid become important and the generalized Stokes relation is no 
longer valid. This scenario is confirmed by detailed calculations using a "volume localization" 
approximation [12]. 

Theory. - In the following, we consider a rotating particle with an oscillating angular 
displacement $(t) ~ $([x))e^*"* where the direction of the vector ^{uj) caracterizes the axis of 
rotation. Since the pressure around the particle stays constant, the oscillating velocity field in a 
pure incompressible Newtonian fluid is described by the Helmholtz equation (V'^+fc^)t)(r, w) = 
with wave number k — ^ iojp/rj [13]. Under stick boundary conditions it assumes the form 

v{r, Lo) = -luj { -) ' [*(^) X r] }—^e'Kr-a) 
\r / 1 — ika 

which determines the external torque Tiuj) = {u})^{ijj) to drive the oscillating particle. 
For small frequencies, the compliance obeys the familiar Stokes result, oi~^ / — iuj = Sira^r]. 
Deviations from this law, occur around ujo = 2r7/(/9a^), i.e., when the penetration depth 
S = Im(fc) equals the particle radius a. 

We now study the equivalent problem for a model viscoelastic medium, represented by a 
two- fluid model [11]: 

= fiV^u+{X + fi)ViV ■u)+T{v-^) (2) 

p-^v = -\7p + T]\7'^v -T{v - —) , div?; = . (3) 

Here an incompressible Newtonian fluid with shear viscosity 77 is coupled to an elastic medium 
with Lame constants A, /i via a conventional friction term. By dimensional analysis, the 
friction coefficient F = r?/^^ contains a characteristic length ^ which is on the order of the 
mesh size in, e.g., an actin network [12]. We neglect the mass density of the elastic medium 
to the one of the Newtonian fluid right from the beginning. The characteristic parameters of 
the theory are the reduced mesh size ^/a, uje = m/?/ and ujq — 2i]/{pa?). The two frequencies 
quantify the fluid inertia, shear elasticity, and shear viscosity. Typical numbers for an actin 
solution [5] are uje — 10^ Hz and wq = 10^ Hz based on a = 3 /^m, -q = O.OIP and ^ = 1 N/m^. 

An oscillating rotating particle creates pure shear fields for both displacement u and 
velocity v (divM = divi; = 0) and keeps the pressure constant. One then derives from Eqs. 

and ^ that the amplitudes u[r,ui) and v(r,w) of the oscillating fields obey a vector 
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Helmholtz equation in analogy to the pure Newtonian fluid, however with replaced by a 
matrix K^. To solve this equation, we introduce 

u{r,uj) = — a^$((jj) X Vg{r) and v{r,Lu) = — (jjea^$(w) x Vh[r) (4) 

and finally arrive at the equivalent Helmholtz equation for g{r) and h{r): 

(V'l + ( ?S ) = with K^ = i( *^/^^, ' \ . (5) 



h{r) J ^2 \^ —iuj/uje 210;^^^ / (ujQa ) — 1 

The reduced radial coordinate is denoted hy f = r/a and the radial part of the Laplace 



operator reads V = :;^^(^^^)- The general solution of Eq. lO is given by exp(il<'r)/r or 



h{r) )~ f \ exp(zyA^f) ^ ^2 / ' ^ ' 

where \i are the eigenvalues of . The matrix S — (ei, 62) is composed of the eigenvectors 
ei and 62 and therefore diagonalizes . Finally, the constants hi are determined by the 
boundary conditions on the surface of the particle. Note that the roots have to be chosen 
such that the exponentials in Eq. ^ decay to zero for large f . 

In general, the external torque on a particle is calculated by T = — J r x crdf where cr is 
the stress tensor and df the directed surface element. In our case, df oc —r and the velocity 
V and displacement vector u point along the azimuthal direction relative to the axis ^(to). 
Therefore, only the respective components of the elastic and viscous stress tensor, 

contribute and give a torque parallel to ^{oj) with magnitude T — T" + T'", where 

= 8^/ia3^H(l + ^[5ii6iAie^^ + 5i262A2e*^]) (8) 
= -iiu8n7]a^^{Lj)(l-l[S2ibiXie'^^ + S22b2X2e'^]) . (9) 



The prefactors on the right-hand side of Eqs. © and © are the Stokes results for a pure 
elastic or viscous medium. 

In the following we consider two cases to determine the constants bi from the boundary 
conditions. In the first case, we assume stick boundary conditions for both the viscous and the 
elastic component at the particle surface, i.e., v = —iuju = —iuj'^{uj) x r\^^^. So we assume 
that the elastic network is attached to the particle. In the second case, the elastic network 
is not attached and the elastic stress tensor cr^^ in Eqs. O vanishes on the particle surface. 
As a result, the elastic torque T" is zero and the displacement field has to fulfill the mixed 
boundary condition {du^p/dr — Utp/r)\^^^ = 0. In both cases we can write down the solutions 
in analytic form. However, the concrete formulas are too large to convey direct information. 
We therefore discuss their graphic representations and determine certain limits. 

Case 1. - In Figs.^andl^we plot the real and imaginary part of the inverse compliance 
a~^{u}), relative to their results for a pure elastic and viscous system, as a function of lu /lo^ and 
Wo/^e! the mesh size is ^/a = 0.1. For low frequencies they both exhibit constant values which 
correspond to the generalized Stokes relation = S'Ka'^Giuj) with G(uj) — fi — iuj-q. This 
result can be extracted from our theory as long as the term 2ui£,-^ / {tuoa^) in the matrix of 
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Fig. 1 - Case 1: Real part of the inverse compliance a ^ , normalized to 9>-Ka^ ^ as a function of uj/(jJe 
and uJo/iOe; the parameter is ^/a = 0.1. 

Fig. 2 - Case 1: Imaginary part of the inverse compliance a~^, normalized to cdSiva^ri as a function 
of Lu/ue and ujo/iUe, the parameter is ^/o = 0.1. 



Eq. (ISJ can be neglected against one, i.e., as long as inertial effects of the fluid are negligible. 
Note that unlike the translational motion, the validity of the Stokes relation extends to w ^ 0; 
there is no lower crossover frequency. The reason is clearly that a rotating sphere only creates 
pure shear fields for both dynamic variables u and v. This indicates that elastic network 
and viscous fluid are strongly coupled to each other and move together. Therefore, their 
dynamics is described by the equation of linear elasticity with the complex shear modulus 
G{oj) (reminiscent to a Voigt element [3]) or alternatively by the Navier-Stokes equation with 
Tj replaced by 77 — ij,/{iuj). With the latter view, the compliance can be calculated using the 
result from the Newtonian fluid [13] with a complex wave number given by k'^ = pw'^ / {fi—iuir}): 



1 + 2Vi + 2x + 2.t3/2/3 - i2x{\ + ^/i)/3 
1 + 2^+ 2x 



1 



Wo 1 + iuJe/oJ 



(10) 



The last equation fits the graphs in Figs.nandlSlwell. It especially accounts for the deviations 
from the Stokes law due to inertial effects at higher frequencies. In a pure Newtonian fluid, 
inertia becomes noticable around the frequency ujq] just set cjg = /i/77 = in Eq. (|10l) . For 
the strongly coupled viscoelastic system, the onset of inertial effects is not so clear from Eq. 
I)10|l . We therefore determined the appropriate crossover frequencies empirically by requiring 
that the compliance a(ijj) deviates from the Stokes law by 10%. As it is already obvious from 
the graphs in Figs. ^ and |21 we flnd that the crossover frequencies exhibit different behavior 
for the real and imaginary part of a~^{uj). For the real part, it scales as ^JuJquJ^ cx p/(pa?) 



whereas for the imaginary part it behaves as ujq 



0.77, ,0.23 



for oJo/uJe < 1 and passes over to loq 



for (jJo/uJe > 1. This cross over is clearly seen in Fig.|21 

So far, we discussed the case of ^/a <C 1. For £,/a ^ 1, the regime of generahzed Stokes law 
still exists but Eq. pO|l does not apply anymore, although the deviations are not dramatic. The 
regime ^/a ~ 1 means that the continuum limit can no longer be applied to elastic networks 
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Fig. 3 - Case 2: Real part of the inverse compliance a ^, normalized to Sira^fi as a function of lu/ujc 
and ^/a; the parameter is luq/cj^ = 100. 

Fig. 4 - Case 2: Imaginary part of the inverse compliance , normalized to Sira'^fi as a function of 
uj/u)e and ^/a; the parameter is wo/we = 100 



such as actin. However, in the two-fluid model ^ just quantifies the frictional couphng between 
fluid and elastic component. One could wonder if there exists a viscoelastic system with such 
a weak coupling so that ^/a ^ 1 is applicable. 

Case 2. - We now address the case where the elastic network is not coupled to the 
particle surface, it only reacts to shear flow via the friction term in Eq. |(5J. In Fig.|3| we plot 
the real part of the inverse compliance as a function of reduced frequency to/uJe and mesh size 
^/a; the additional parameter is set to wo/we = 100. For small frequencies, the real part is 
close to zero, in contrast to case 1, where it assumes the reference Stokes value of Sna^^, as 
already discussed. The friction between the two components is sufficiently small so the fluid 
permeates the elastic network without deforming it noticeably. Then for increasing frequency, 
an edge occurs and Ke{a~^) enters the region where the Stokes law is valid. Note, however, 
that this region is considerably reduced compared to case 1. Correspondingly, the imaginary 
part in Fig. ^ (also in units of Sna^^l) exhibits a ridge. This is what we roughly expect since 
real and imaginary part of the compliance are connected via Kramers-Kronig relations [14]. 
The features described so far can be explained by the formula 



U^T^ { 1 

icjT ( 1 

1 — iu^T \ UJ„T 



1 a 



with T = — — v^l - iuj/Uf, , (11) 
We oc; 



which follows from our theory when we neglect again the "inertial term" 2^0^^"^ / {ujqo?) in the 
matrix of Eq. (0) and set ^/a <C 1. To determine r, the complex root with positive real 
part has to be taken. For lo/uo,, <C 1, the existence of the edge and ridge in relation Hll(l 
are obvious. Furthermore, the relation demonstrates that edge and ridge, or the onset of the 
Stokes regime, occur at a frequency which scales as ^/a. This is in contrast to translational 
motion where it scales as (C/a)^ [5, 12]. In Fig.Elthe imaginary part of with proper scale 
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Fig. 5 - Case 2: Imaginary part of the inverse compliance q ^, normalized to a;87ra^r; as a function 
of LojijJe, and ^/a; the parameter is ujo/uje ~ 100 

factor uSna^ri is plotted. It reveals the strong effect of friction between both components 
when the fluid permeates the elastic network. Note the logarithmic scale of the vertical axis. 
Inertial effects for ^/a <C 1 are again described by Eq. (^U)). So the Stokes regime in Fig. |5l 
extends up to uj/uie = 10 in agreement with Fig. El Nevertheless, we find that in case 2, the 
validity of the Stokes law is reduced. 

Conclusions. - Based on a two-fluid model, we determined and discussed the compliance 
for a rotating spherical particle embedded in a viscoelastic medium. In the case where both 
the viscous and elastic component obey stick boundary conditions at the particle surface, we 
identify the validity of a generalized Stokes law, and therefore a simple relation to the complex 
shear modulus, starting from zero frequency and limited only at high frequencies by inertial 
effects. This is in contrast to translational motion. When the elastic network is not coupled 
to the particle, the compliance exhibits a region of low elasticity and high effective viscosity 
starting from zero frequency indicating the strong friction which occurs when fluid permeates 
the elastic network. The regime of the Stokes law is reduced. 

Real viscoelastic systems such as actin solutions possess a frequency-depended complex 
shear modulus [2,4,5]. So to discuss the frequency dependence of the compliance, e.g., in 
case 1, one has to consider a non-trivial path in Figs. ^ and [21 determined by the varying 
parameter uJo/uje oc rf' (uj) / . This is, however, no problem in the Stokes regime but leads 
to additional effects outside the Stokes regime as, e.g., in case 2. 

Our study of different boundary conditions clarifies that the interpretation of compliances 
measured in experiment needs care. So far, in the translational case, always stick-boundary 
conditions are assumed. Our results show that some slip of the elastic network changes the 
measurable compliance dramatically. This could lead to false interpretations. 

Clearly, one-bead microrheology with rotating particles also suffers the drawback that it 
does not probe bulk properties, as it is done with two-bead methods [15]. On the other hand, 
since in the Stokes regime the shear fields decay as 1/r^ instead of 1 /r for translational motion, 
the method could be used to monitor explicitely the effect of the embedded particle on the 
elastic network [16]. A challenge will be to develop the theory for two-bead microrheology 
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including the rotational degree of freedom [17]. 

In this article we laid the theoretical basis for microrheology with rotating particles. Our 
results demonstrate that it can be a useful extension of existing methods based on translational 
motion and also complements them. So we hope that our work stimulates further experimen- 
tal investigations which use rotating particles as colloidal probes to explore the mechanical 
properties of soft materials, especially in connection with biological systems. 
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